Entropy of gravitationally collapsing matter in FRW universe models 
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We look at a gas of dust and investigate how its entropy evolves with time under a spherically 
symmetric gravitational collapse. We treat the problem perturbatively and find that the classical 
thermodynamic entropy does actually increase to first order when one allows for gravitational po- 
tential energy to be transferred to thermal energy during the collapse. Thus, in this situation there 
is no need to resort to the introduction of an intrinsic gravitational entropy in order to satisfy the 
second law of thermodynamics. 
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I. INTRODUCTION 

Today there is broad consensus among cosmologists 
that the configuration of energy in the early universe 
was very homogeneous and isotropic. Observations of 
the temperature variations in the cosmic microwave ra- 
diation have also shown that the Universe was in a state 
close to thermodynamic equilibrium 400 000 years after 
Big Bang, with relative temperature and density varia- 
tions of the order 10 -5 0. 

Naively, one expects the entropy in a gas to be higher 
the more homogeneously distributed its density and tem- 
perature is. Thus, the early universe described above 
should be one of near maximal entropy, since it differs 
only by a small fraction from one of total homogeneity 
in density and temperature. However, due to gravity, 
small inhomogeneities start to grow and eventually end 
up forming structures such as galaxies, stars, planets, 
planetary clouds etc. This evolution is in the direction of 
greater inhomogeneities both in energy density and tem- 
perature, which according to the argument above, ap- 
pears to violate the second law of thermodynamics by 
decreasing the entropy. Obviously, something must be 
wrong with this picture, since we consider the second 
law of thermodynamics to be a basic law of physics and 
it should therefore not be violated. 

A possible solution to this apparent paradox comes 
from considering the quantity known as gravitational en- 
tropy. This was introduced by R. Penrose in the 1977 in 
connection with his study of the properties of the initial 
singularity of the universe It is a quantity which 

can be interpreted as an entropy intrinsic to the gravita- 
tional field. It takes into account the attractive nature 
of gravity and increases as a gas collapses under the in- 
fluence of gravity. This allows one to define a general 
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entropy which is the sum of the ordinary thermodynamic 
entropy and this new gravitational entropy. If the sum 
of the two types of entropy increases during gravitational 
collapse of a gas, the second law of thermodynamics will 
then be preserved. 

In this paper we will show that the thermodynamics 
entropy of a collapsing gas does actually increase, which 
allows us to explain the collapse without introducing the 
gravitational entropy. We look at a perturbed ideal gas in 
a FRW background and consider changes in its classical 
entropy up to first order in the energy density. We find 
that the increase in the thermal energy which comes from 
potential energy released in the collapse actually makes 
the total thermodynamic entropy increase, even though 
the temperature inhomogeneity increases. 

The structure of this paper is as follows. In section 2 
we look at a simplified model consisting of ideal particles 
in a box and explain why one would expect the ther- 
modynamic entropy to decrease as the inhomogeneities 
increase. In section 3 we introduce a tool which we 
will need when considering the growth of small inhomo- 
geneities, namely cosmological perturbation theory. In 
section 4 we derive an expression for the thermodynamic 
entropy of a gas in an expanding universe. In section 
5 we specialize to spherically symmetric collapsing gases 
and arrive at our main result. Finally, section 6 contains 
a summary and our conclusion. 



II. SIMPLE PICTURE: IDEAL GAS IN A BOX 

In this section we look at gas confined to a box and 
show that its entropy is maximal when the density is ho- 
mogeneous and the temperature is the same everywhere. 

Consider an isolated box that is divided into two cham- 
bers of equal volume. Each of these contains a gas of 
the same type of particles with different temperatures 
and densities, as illustrated in Fig. ^ We will look at 
two different scenarios: 1) when the temperature in the 
chambers is the same but the density is different. And 
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Figure 1: An isolated box consisting of two separated cham- 
bers with particles of different temperatures and densities. 



2) when temperature is different, while the density is the 
same. We then compare these to a third scenario, in 
which we remove the wall between the two chambers so 
that we have only one gas with just one temperature and 
one density. 

First, we need the expression for the entropy of an ideal 
gas consisting of N ideal particles in a volume V 
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Let us look at the third scenario first. Let T be the 
temperature in the gas and N the total amount of parti- 
cles in the box. The entropy in the box in this scenario, 
which represents the totally homogeneous case, is then 



S 3 = Nk B In 
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where we have defined the constant 

r . "mfc B Te 5 / 3l3/2 
K = 



2nh 2 



(2) 



(3) 



How does the entropy differ from this in the other two, 
inhomogeneous scenarios? In the first scenario the tem- 
peratures are the same, but the number of particles in 
the two chambers is different. The total number of par- 
ticles is conserved, so we have that Ni + N 2 = N, where 
N\ and N 2 are the numbers of particles in the two cham- 
bers respectively. Furthermore, we assume that the two 
chambers have the same volume V/2. The total entropy 
in the box is then the sum of the entropies in the two 
chambers: 
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A transition from the totally homogeneous scenario to 
this results in an entropy difference which is 



AS = Si- S 3 = -Nk B \^ln2+x lnx+(l-x) In (1 - x) 

(5) 

where we have defined x = Ni/N. We can restrict x to 
the interval < x < 1/2 without any loss of generality. 
It is then a simple task to verify that AS < 0. Thus, 
based on this simple picture, we can conclude that an 
increase in inhomogeneity of the density of an ideal gas 
leads to a reduction of the entropy. 

Now, let us look at the second scenario, where there 
is a temperature difference between the two chambers 
which both contain the same amount of particles, N/2. 
The temperatures of the two chambers are T\ and T 2 
respectively. The total entropy of the box in this scenario 
is 



S 2 = Nkt 
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If we imagine the gas first being in the totally homo- 
geneous state of scenario 3 and then changing into the 
thermally inhomogeneous state of scenario 2, the entropy 
difference will be 



AS" 



3Affc S] T X T 2 
In r^— 

4 T 2 



(7) 



Assume now that the thermal energy is conserved in this 
transition, i.e. that the average temperatures in the two 
scenarios are the same. This means that 7\ + T 2 = 2T. 
Under this assumption, when does this entropy difference 
become non-negative? We see that this will be the case 
when 



4TiT 2 



(Ti + T 2 f 



>1 i.e. (Ti-T 2 ) 2 <0. 



(8) 



This is only satisfied when T\ — T 2 = T, in which case 
there will be an equality between the left and the right 
hand side. If Ti and T 2 are different, i.e. there is a tem- 
perature difference between the two chambers, the en- 
tropy difference in JJJ will be negative. Thus, an increase 
in temperature inhomogeneity will result in a decrease in 
the entropy. 

However, we must not forget that we have assumed 
that the thermal energy is conserved, just as we assumed 
that the particle number is conserved. For this simple ex- 
ample of particles in an isolated box both these assump- 
tions will be true. For a gravitationally collapsing gas, 
however, this need not be true. The total mass, which is 
the equivalent of the total particle number in the box ex- 
ample, must obviously be conserved so we should expect 
the entropy to decrease as the energy density becomes 
more inhomogeneous. But the thermal energy will not be 
conserved as the gas undergoes a gravitational collapse. 
As the energy density near the overdensity increases, the 
potential energy of the inward falling portions of the gas 
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is converted to thermal energy. In other words, there will 
be an increase in the thermal energy of the gas, which 
tends to increase the entropy. 

To summarize, the entropy change in a gravitationally 
collapsing gas can be ascribed to two different effects, 
namely a decrease due to increasing density and temper- 
ature inhomogeneities and an increase due to increasing 
temperature. As we will show, using first order pertur- 
bation theory, the sum of these two effects yields an in- 
creasing total entropy when one assumes that all the loss 
in potential energy is converted into thermal energy. 

We start by reviewing the basics of scalar perturbation 
theory. 



III. SCALAR PERTURBATION THEORY 

We will only consider scalar perturbations since these 
are the only ones that give rise to gravitational collapse. 
For a more detailed review of perturbation theory, the 
interested reader is referred to the standard references 
[E IE @- We assume that the universe is occupied by 
matter in the form of a perfect fluid with no anisotropic 
stress. The energy density inhomogeneity is described as 
a linear perturbation to a flat, matter dominated FRW 
universe. This allows us to write the perturbed metric 
in terms of only one perturbing function $, the so-called 
Bardeen potential: 



ds 2 



\r)) {(1 + 2$)d?7 2 - (I - 2<^)5 i] dx l dx 1 } , (9) 



where 77 is conformal time and a is the scale factor of the 
universe. 

The energy-momentum tensor for the matter content 
is written as a homogeneous zeroth order term plus a 
non-homogeneous first order perturbation: 



6T£. 



(10) 



Using the definition for the energy-momentum tensor of 
a perfect fluid, the equation of state for matter, and the 
four- velocity identity u^u^ = I, we can write the compo- 
nents as: 
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where po is the average density of the fluid, Sp is the 
density perturbation and Su 1 the velocity perturbation. 
In the expressions above and throughout this paper we 
use the convention that Latin indices run over spatial 
components only, while Greek indices run over both space 
and time components. 

We require the Einstein equation for the fluid to be 
satisfied independently for each order in the perturba- 
tion. This gives us the following zeroth order equations 



H 2_ 8 „^2. 
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(14) 



and 

H 2 + 2H = 0, (15) 

where H = i jj 2 and the dot denotes differentiation with 
respect to rj. The first order equations are 

3, 



V 2 $ - 3H($ + = -H 2 S 
\$ + H®} = -^H 2 aSu % 



(16) 
(17) 
(18) 



$ + 3W$ + (H 2 + 2H)§ = 

where we have defined the density contrast as 8 = 5p/po, 
and A t j = J!4-. The zeroth order equations are the or- 
dinary FRW equations for a matter dominated, flat uni- 
verse expressed in conformal time instead of the usual 
comoving time. The solution to these equations is 



a = I — ) and p = Poo [ — 



V 



(19) 



where we have defined rjo such that a (770) = 1, and poo is 
the energy density at rj — 770, which can be written as 



Poo 



(20) 
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We shall later need the relation between conformal and 
comoving time. Using the definition a{rj)drj — dt, we can 
write this as 



(21) 



where rj is the initial comoving time that corresponds to 
r/Q. These are related via the expression 

m 



(22) 



Let us introduce a new dimensionless time parameter r 
which measures time relative to the initial time to, i.e. 

t - 1 



to 



(23) 



Using this new time parameter we can write the scale 
factor and the unperturbed energy density as 



o(r) = (l + r)3 and p (r) = p o(l + t)~ 



(24) 



Next, we solve the first order equations. We start with 
equation l|18|) . which doesn't couple to the other two 
equations, and obtain the metric perturbation <i>. The 
remaining perturbing functions S and 5u l are then ob- 
tained by a simple substitution of $ into equations i|16|) 
and (|17|) . Disregarding solutions that decrease with time, 
we can write the perturbations as: 



$(x,t)=/(x) 



(25) 



*(x, r) = 1 % 2 (1 + r) 2 / 3 V 2 /(x) - 2/(x) (26) 



^(x,r) = -|^(l + r)- 1 /» 



(27) 
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where /(x) is an arbitrary function of spatial coordinates 
only. In order for these perturbations to be physically 
acceptable they must vanish at infinity, i.e. 



lim /(x) = lim 

x|| — >oc ||x|| — >o 



5(x, t) = lim <5u J (x, r) = . 



(28) 

Furthermore, we must also require that the total energy 
density at t = in the perturbed and the unperturbed 
universe remains the same, which is the same as saying 
that the total energy must be conserved when the pertur- 
bation is introduced. For this requirement to be satisfied, 
the initial density perturbation must satisfy the integral 
condition 



/ 5{x,0) dV = 0. 
Jv 



(29) 



Substituting for the left hand side from Eq. 1261) and 
using the boundary conditions in Eq. I|28|) . we find that 
the volume integral of the metric perturbation must also 
vanish, 



f f(x)dV= f V 2 f(x)dV 
Jv Jv 



0. 



(30) 



This implies that the volume integral of the density per- 
turbation must vanish for all values of t, 



/ 5(x,r) dV 
Jv 



(31) 



What this equation says is simply that the total energy 
in the perturbed universe must be the same as that in 
the unperturbed universe at all times. This is nothing 
but a statement of energy conservation. 



IV. ENTROPY OF A PERTURBED IDEAL GAS 
IN A FRW UNIVERSE 

Equation gives us the entropy of an ideal gas con- 
sisting of N distinct particles. The collapsing gas we 
wish to examine, whose time evolution is given by Eq. 
(|26J) . consists of a continuous fluid. We must therefore 
rewrite the expression in Q into a form that we can use 
for a continuous fluid. In order to do that we consider an 
ideal gas contained within a small volume element dV. 
The number of particles inside this volume is 



dN = 



pdV 



(32) 



where m is the mass of the particles which the gas consists 
of. Inserting this expression for the particle number into 
Eq. 0, we can write the entropy associated with the 
volume element in terms of the density of the fluid: 



dS = kn — In 



dV = adV . 



(33) 



where the constant K is defined in J3J) and a can be 
interpreted as the entropy density of the ideal, continuous 
gas distribution. 

We substitute the energy density of the perturbed pres- 
sureless gas for the density which appears in this expres- 
sion. The former can be written as p — /Oo(l + S), where 
po and 8 are given by l|24|) and (|26J) respectively. This 
allows us to write the entropy density as 



KT 3/2 



PO 
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, mKT 3 / 2 
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Time dependence enters into this expression via the un- 
perturbed energy density po, the density contrast 5 and 
the temperature T. For a totally homogeneous uni- 
verse which contains only matter, the temperature can 
be shown 9] to scale like T ~ a~ 2 , where the bar de- 
notes that the temperature is that of a non-perturbed 
gas. In terms of the dimensionless time parameter r, 
we can write the time dependence of the homogeneous 
temperature as 



T = T (l + r) 



-4/3 



(35) 



where To is the temperature of the gas at the initial time 
t = 0. 

In a perturbed gas we expect there to be an addi- 
tional, non-homogeneous contribution to this tempera- 
ture. Thus, we can write the total temperature as 



T = f(l + A r ), 



(36) 



where TAt is the non-homogeneous addition to the ho- 
mogeneous temperature T resulting from the first order 
density perturbation 5. As we will see later, At will 
turn out be too large for us to treat it as a first order 
perturbation. 

The time evolution of At depends on how much en- 
ergy we assume is transferred from potential into thermal 
energy due to the gravitational collapse and how this 
is transferred. The two extremes are: 1) No energy is 
transferred and 2) All the potential energy is transferred 
adiabatically into thermal energy. The most realistic sce- 
nario would probably be somewhere in between these two 
extremes, but for simplicity we will assume the latter sce- 
nario when we calculate At explicitly in the next section. 
Using Eqs. and (|31)|) we can write the entropy 

density of the ideal, perturbed gas as 
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The total entropy of the gas is the volume integral of this 
expression over the whole of space. The volume element 
which appears in the integral is given by the determinant 
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of the spatial metric hij . 
can be written as 



In Cartesian coordinates, this 



dV = ^\deth l0 \d 3 x = (1 + r) 2 (l - 3<S>)d 3 x , (38) 

where d 3 x is the Euclidean volume element. This allows 
us to write the entropy element as 
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kBPoo 
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(39) 



As a consistency check we can calculate the entropy 
inside a comoving volume V of the unperturbed FRW 
model, which we know to be a constant Using our 
definition of the entropy of an ideal, cosmological gas 
ijnij , we find that 
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which is indeed a constant. What about the perturbed 
entropy? If the gravitational collapse is not to conflict 
with the second law of thermodynamics this should in- 
crease with time, or at least not decrease with time. The 
change in entropy resulting from the density perturbation 
is 



AS = 



kBPoo 



<J(ln 
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d 3 x. 
(41) 



Due to the energy conservation equations H3()[l and 1)3 l|l. 
the first two terms in this integral will vanish. This leaves 
us with a contribution from the temperature inhomo- 
geneity only 



As = i k B poo 
2 m 
_ 3 k B poo 
~ 2 m 



f (1-3$) ln(l + A T )d 3 
Jv 

In (1 + A T ) d 3 x . 



(42) 



In the simplified box example of section[n]we found that 
density inhomogeneities tend to reduce the total entropy. 
It might therefore seem odd that the density inhomo- 
geneities in the gravitationally collapsing gas don't con- 
tribute to the entropy. The reason that the entropy de- 
creased in the box example is that the expression for the 
entropy is non-linear in the particle number (or alterna- 
tively in the density). This leads to the result that the 



sum of the entropies of two gases with different particle 
numbers is generally different from that of a gas whose 
particle number is equal to the sum of the particles in 
the two first gases. But if we treat the problem pertur- 
batively to first order, we force all physical quantities to 
be linear in the perturbed quantities. Since the density is 
conserved, as we can see in Eq. 1)3 ljl . the contribution to 
the entropy from the density inhomogeneity must vanish 
to first order. 

In the next section we examine how the entropy 
changes when energy is converted adiabatically from 
gravitational potential energy into thermal energy for a 
spherically symmetric collapsing gas. 



V. SPHERICALLY SYMMETRIC 
PERTURBATIONS 

In this section we will specialize to spherically symmet- 
ric perturbations. Furthermore we assume that the po- 
tential energy in the gravitational field is transferred into 
thermal energy adiabatically during the collapse. The 
density inhomogeneity will then give rise to a tempera- 
ture gradient. If we assume that the gas is in hydrostatic 
equilibrium during the whole collapse and that the gas 
collapses towards the origin of the coordinate system, we 
can write the temperature gradient as |ll| 



d At ^AirGm pa 2 
dr ks Tr 2 



dr'r' 2 5(r',T). 



(43) 



where C = 2/5 for an ideal monoatomic gas. To proceed 
further we need to know how the density perturbation <5 
behaves. We must therefore solve the differential equa- 
tions <|25|l and (|26|l . In order to do this we must specify 
initial and boundary conditions for both the metric and 
the density perturbation. The boundary conditions are 
stated in Eq. (J2SJ. Once we give an initial profile for the 
density perturbation we can solve Eq. at r = for 
the metric perturbation. The time evolution of 5 is then 
determined by reinserting the solution for $ into i)26[) for 
arbitrary r. 

We restrict ourselves to a general type of initial condi- 
tions where there is only one initial overdensity, namely 
centered at the origin. These profiles must satisfy the 
energy conservation condition in l|29|) . A simple density 
profile contained within this class of initial conditions is 



S(r, 0) = d 



r 
L 



1 /rx = 
3 \L 



(44) 



where L and do are measures of the size and the ampli- 
tude, respectively, of the initial overdensity. The reason 
that we have chosen this explicit expression for the initial 
density perturbation is that it allows us to solve the dif- 
ferential equation analytically. However, as we will show 
in the appendix, the results we obtain apply qualitatively 
for all initial density profiles in the class we defined above. 
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In analogy with the time parameter r, we introduce a 
new dimensionless radial coordinate y, which measures 
comoving radial distance relative to the length scale L, 



y =L 



(45) 



The differential equation which determines the metric 
perturbation can now be written as 



2 d 
V dy 



do 



1 



y 



(46) 



This can be solved analytically by e.g. using the com- 
puter program Maple. However, the analytical solution 
is too long for us to list up here. Instead we illustrate 



the solution by plotting it in Fig. 2(a) The amplitude 
of the initial perturbation was chosen to be do = 10~ 5 . 
This corresponds to the amplitude of the density pertur- 
bations in our own universe at the time of recombination, 
to = 400 000, which we choose as the initial time of our 
perturbation. 

The time evolution of the density perturbation de- 
pends on the size of the perturbation. If L is sufficiently 
large, the ratio ^ will be so small that S(y, r) remains 
essentially constant in time, which can seen directly from 
Eq. (|2oTl . Since we are interested in perturbations that 
grow with time, L must be chosen accordingly. This can 

> 1. The value we used to 



be achieved by choosing ^ 
obtain the solution plotted in Fig. 2(a) 



was 



10. 



According to Eq. (|26(l . once we know the metric per- 
turbation we automatically know the time evolution of 
the density perturbation. The analytical expression for 
this is even longer than that for the metric perturbation 
so we will omit writing it down. Fig. 2(b) shows a plot 
of the density contrast for a selection of different times, 
thus illustrating how it grows with time. 

Obtaining the relative temperature change At is now 
simply a matter of integrating Eq. i|43|) with the density 
contrast given by the analytical expression found above. 
Just as for the two perturbed quantities 6 and <!>, we 
demand that the At vanishes as r — > oo. This allows us 
to write the relative temperature change as 



A T (j/,r) = - 



6mC /LV/, ^ 
— — - (l+r)> 
fcs Too \W 



/V4> f dy"y" 2 5(y",r). 

Joe y Jo 



(47) 



The remaining constant which we need to determine in 
this expression is the initial temperature Too- Since the 
initial time and amplitude of the density perturbation 
were chosen to correspond to perturbations in our own 
universe at the time of recombination, it is only natu- 
ral that we also choose T 00 to be the temperature of the 
universe at the same time, namely Too = 3 000 K. Fur- 
thermore, the gas we consider consists of baryons, which 




(a) 




■ 1=0.0 

■ 1=0.25 
1=0.5 
1=0.75 
1=1.0 



(b) 



Figure 2: (a) The left panel shows a plot of the metric per- 
turbation that results from the initial density inhomogene- 
ity described in Eq. 1441 . (b) The right panel illustrates 
the time evolution of the corresponding density perturbation. 
The curves represent the spatial configuration of the density 
contrast at, from bottom to top, r = 0.0, 0.25, 0.5, 0.75, 1.0 



allows us to use the mass m = 1.67 x 10~ 27 kg. Us- 
ing these values we find that the dimensionless constant 
that multiplies the integrals in Eq. I|47|l will be of the 
order ~ 10 7 . Since the density contrast is of the order 
~ 10~ 5 we see that the relative change in temperature 
caused by the density perturbation must be very large. 
We can find a plot of this relative temperature change 
in Fig |21 which shows us that At is positive everywhere 
and that it grows with time. Thus, inserting the ana- 
lytical expression for At into Eq. (|42|l . we see that the 
entropy change induced by the density perturbation will 
be positive, and it will grow with time as the density in- 
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Figure 3: A plot of the temperature perturbation th at cor rc- 
sponds to the density perturbation illustrated in Fig|2(b)| 



homogeneity increases. This shows that the entropy of 
a gravitationally collapsing gas evolves according to the 
second law of thermodynamics up to first order, ft is 
natural to assume that this will be the case to any order. 

Strictly speaking, we have only showed this for the 
special initial density perturbation described in Eq. 
fn the appendix we show that this will be the case for all 
initial density perturbations of the class defined in the 
beginning of this section. 



VI. CONCLUSIONS 



The main objective of this paper was to show that the 
entropy of an inhomogeneous gas increases as the gas col- 
lapses under the influence of gravity. Naively, one might 
expect the opposite to be the case since inhomogeneities 
in both the density and the temperature increase under 
such a collapse, which is an evolution that we generally 
associate with a decrease in entropy. By allowing for 
a transfer of energy from the gravitational potential to 
thermal energy, the temperature in the gas will increase 
as a result of the collapse. Treating the inhomogeneous 
gas as a first order perturbation to a homogeneous FRW 
model, we showed that the increase in temperature re- 
sults in an increase in the entropy which outweighs any 
decrease due to increasing inhomogeneities in the tem- 
perature and and the density. This was shown to be the 
case for any initial density inhomogeneity which consists 
of one overdense region. Although our results were de- 
rived only up to first order in the inhomogeneity, it is 
only natural to extend the conclusions to any type of 
inhomogeneity. This allows us to conjecture that the en- 
tropy of a any gravitationally collapsing gas will always 
increase with time, in accordance with the second law of 
thermodynamics. 



We wish to thank professor K.P. Tod for a very useful 
correpsondence 



Appendix: MORE GENERAL PERTURBATIONS 

We define a class of perturbations where the initial 
profile has only one overdensity. The coordinate system 
is chosen so that the center of the overdensity is situated 
at the origin. Define the function 



F(y,r) 



From Eq. (|ol"|l we get that 



dy'y l2 5(y',T) 



lim F(y, r) = . 

y— >oo 



(A.l) 



(A.2) 



The fact that the only overdensity is that at the origin 
implies that there exists a yo such that S(y, r) > for 
< y < 2/o and d(y,r) < for y > yo- This along with 
Eq. HA.2p implies that F(y, r) must be positive for all 
values of y, 



F(y, t) > for y > . 



(A.3) 



The temperature change in Eq. I|47[l can be written as 



A T (y,r) oc -(1 + r) 



2/3 



y 2 



(A.4) 



From this expression we see why the lower limit of the 
integration must indeed be +00: We require the temper- 
ature change to vanish at infinity. Since the integrand is 
always greater than zero, this can only be accomplished 
if the range of integration vanishes at infinity, which im- 
plies that the lower integration limit must be +00. This 
in turn implies that the integral in Eq. (|A.4|I will be 
negative for all y. Thus we have shown that At will be 
positive for all y. This is the first step in our proof. We 
must also show that At increases with time. In order to 
do this we differentiate Eq. <|A.4|) with respect to r. This 
yields the expression 



dAr 



oc 



+ (l + r) 2 / 3 



y dy' dF(y' 7 r) 



y 



0t 



(A.5) 



The first term inside the parentheses will be negative due 
to the same arguments as above. Using the definition in 
Eq. (H3J, we can write 



H(y,r) 



_ dF(y,r) 
dr 



= jdy>y" d -^l. (A.6) 



Again, using Eq. I|A.2|) . we find that lim^^oo H(y, r) = 0. 
We know that the effect of gravity on the density pertur- 
bation is such that overdense regions become more dense, 
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while underdense regions become less dense. This means 
that > for < y < y and < for y > y . 

Thus, just as for the integrand in Eq. 1A.4|) . this implies 
that the integrand in Eq. i|A.6|) must be positive and 
hence that H(y, r) must also be positive. The integral in 
the second term in Eq. ljA.51) must therefore be negative 
since the integrand is positive while the integration path 



is negative. This proves that is positive for all y > 
and r > 0. 

In summary, we have shown that a density perturba- 
tion of the class defined at the start of this appendix 
yields a temperature change At which is positive every- 
where and grows with time. This concludes our proof. 
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